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ABSTRACT 

We report on the detection of four rotating massive disks in two regions of 
high-mass star formation. The disks are perpendicular to known bipolar outflows 
and turn out to be unstable but long lived. We infer that accretion onto the em- 
bedded (proto)stars must proceed through the disks with rates of ~10~^ yr~^. 

Subject headings: Stars: formation - Radio lines: ISM - ISM: molecules - ISM: 
individual (G24.78+0.08, G31.41+0.31) 



1. Introduction 

The formation of massive stars represents a puzzle from a theoretical point of view. 
Unlike their low-mass counterparts, they are believed to reach the zero-age main sequence 
still deeply embedded in their parental cores: in particular, Palla & Stahler (1993) predict 
that this occurs for stellar masses in excess of 8 Mq. Once the star has ignited hydrogen 
burning, further accretion should be inhibited by radiation pressure and powerful stellar 
winds, with the consequence that stars more massive than 8 Mq should not exist. Two 
formation scenarios have been proposed to solve this paradox: non-spherical accretion (Yorke 
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& Sonnhalter 2002) and merging of lower mass stars (Bonnell, Bate, & Zinnecker 1998). 
Discriminating between these two models represents a challenging observational goal. 

In this context, the detection of disks would strongly favour the accretion scenario, since 
random encounters between merging stars are not expected to lead to axially symmetric 
structures. On the contrary, conservation of angular momentum is bound to cause flattening 
and rotation of the infalling material, thus producing disk-like bodies. Indeed, circumstellar 
disks have been detected in low-mass stars and found to undergo Keplerian rotation (Simon, 
Dutrey, & Guilloteau 2000). Similar evidence has been found in a few high-mass young 
stellar objects (YSOs), but in most cases the angular resohition was insufficient to assess 
the presence of a disk unambiguously. In conclusion, only few bona fide examples are known 
(see Cesaroni 2002) and all of these are associated with moderately massive stars (Bl-BO). 
This is not sufficient to understand the role of disks in the formation of even more massive 
stars and estabhsh the relevance of accretion to this process. 

With this in mind, we have decided to perform a search for disks in a limited number 
of high-mass YSOs. For this purpose, we have selected two luminous objects with typical 
signposts of massive star formation such as water masers and ultracompact (UC) Hll regions. 
The first, G31. 41-1-0. 31 (hereafter G31), is a well studied hot core located at 7.9 kpc (Olmi 
et al. 1996b; Cesaroni et al. 1998), where preliminary evidence of a rotating massive disk 
oriented perpendicularly to a bipolar outflow has been reported in Cesaroni et al. (1994). 
The other, C24.78+0.08 (hereafter G24), is a cluster of massive (proto)stars with a distance 
of 7.7 kpc, where recently Furuya et al. (2002) have detected a pair of cores, each of these 
associated with a compact bipolar outflow. By analogy with G31.41 the expectation is that 
also in this case the cores could contain rotating disks perpendicular to the flow axes. 

On the basis of previous experience with this type of objects (see, e.g., Cesaroni et al. 
1999), CH3CN has been used as disk tracer. This is a low-abundance molecule which is 
excited in very dense regions. Therefore, searching for disks requires not only high angular 
resolution, but also great sensitivity given the faintness of the fines observed. In order to 
achieve these goals, we have used the Plateau de Bure interferometer (PdBI) at 1.4 mm in 
the most extended configuration. In this letter we present the discovery of four rotating disks 
in the two regions studied and discuss possible implications for the star formation process. A 
full report on the results obtained and a more detailed analysis of the data will be presented 
in a forthcoming article. 
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2. 



Observations 



We carried out observations in the 1.4 mm continuum and CH3CN(12-11) line emission 
with the PdBI on 2003 March 16. The inner hole in the {u,v) plane has a radius of 15 kX. 
Line data have been smoothed to a spectral resolution of 0.5 km s^^ and channel maps were 
created with natural weighting, attaining a resolution of l''2x0'.'5 (full width at half power 
of the synthesized beam) and a sensitivity of 40-50 mJy /beam/channel {la RMS). For the 
continuum map the resolution and sensitivity are l'.'2x0'/5 and 4-6 mJy/beam, respectively. 

In the following, we further analyze previous CS(3-2) observations obtained with the 
Nobeyama Millimeter Array (NMA) by Cesaroni et al. (2003), to whom we refer for technical 
details. 



The main goal of our study was the discovery of rotating disks associated with massive 
YSOs deeply embedded in dense, compact cores. This was achieved searching for well defined 
velocity gradients in the cores, perpendicular to molecular outflows powered by the YSOs, as 
illustrated in the next sections. 



The region G31 consists of a hot core, detected in various high-energy lines (Olmi, 
Cesaroni, & Walmsley 1996a) located at the center of a bipolar outfiow, at ~5" from an UC 
Hll region (Cesaroni et al. 1998). On the other hand, G24 is more complex, as it contains 
four distinct objects (see Fig. 1 of Furuya et al. 2002): two of these, G24 A and G24 C, are 
massive cores associated with two bipolar outfiows and represent the target of the present 
study. 

A picture of the G24 A and G31 cores is given in Fig. 1, where overlays of the 1.4 mm 
continuum and integrated CH3CN(12-11) line emission are shown. Note that no map is 
shown for G24 C because no CH3CN(12-11) line emission has been detected with the PdBI. 
However, the 1.4 mm continuum flux is consistent with the extrapolation of the spectral 
energy distribution presented by Furuya et al. (2002), thus confirming the existence of such 
core. Very likely the fact that G24 C is detected in the CH3CN(8-7) transitions (Furuya et 
al. 2002), but not in the (12-11) is due to this core being significantly colder than G24 A 
(see Codella et al. 1997), which makes it difficult to detect high energy fines. 



3. 



Results and discussion 



3.1. 



Structure of the cores 
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When observed with sub-arcsec resolution, G24 A is resolved into two separate cores. 
This is evident both in the 1.4 mm continuum and line maps. In the following we shall refer 
to these cores as G24 Al (the one to the SE) and G24 A2 (the one to the NW). The former 
lies slightly closer to the geometrical center of the bipolar outflow reported by Furuya et 
al. (2002), but the small separation between the cores and the fact that they arc aligned 
along the outflow axis make it difficult to establish whether the outflow is indeed associated 
with G24 Al; in the following we arbitrarily assume that this is the case. The conclusions 
derived in our study are independent of the association of the outflow with either of the 
cores. Noticeably, G24 Al coincides with an unresolved UG Hll region detected by Godella, 
Testi, & Gesaroni (1997), whereas no free- free emission is reported towards G24 A2. To 
determine whether this is an effect of different evolutionary stages of the two cores requires a 
detailed comparison of their physical properties which we postpone to a forthcoming paper. 

The appearance of G31 is even more intriguing: while the 1.4 mm continuum seems 
to trace a roughly spherical core, the GH3GN map reveals a toroidal structure with the dip 
centered at the position of the continuum peak (Fig. 1). This suggests that cither the GH3GN 
abundance drops dramatically in the central region of the core, due to a temperature increase 
towards the center (in agreement with the findings of Olmi et al. 1996b), or the temperature 
near to the embedded source is so high that the ground level states of CH3CN are under 
populated. Interestingly, the two peaks of the GH3GN emission are roughly symmetric with 
respect to the axis of the outflow observed by Olmi et al. (1996b), suggesting a physical 
connection between the toroid and the flow. Such a connection will become more evident 
when considering the velocity field in the cores. 



3.2. Kinematics of the cores 

An obvious way to analyze the velocity field in the cores is to produce maps of the 
line peak velocity obtained with Gaussian fits. Since multiple GH3GN X-components are 
simultaneously observed in the same intermediate-frequency bandwidth, it is possible to 
improve the accuracy of the fit by fitting all lines together, assuming identical widths and 
fixing their separations to the laboratory values (see, e.g., Olmi et al. 1993). Such a fit 
has been made in each point where CH3CN emission is detected. The maps of the LSR 
velocity for the G24 Al, G24 A2, and G31 cores can be seen in Figs. 2c, 2d, and 2e. The 
same method could not be applied to G24 C because the CH3CN(12-11) line emission is not 
detected towards this core in our PdBI observations, while the spectral resolution used for 
the GH3GN(8-7) transition by Furuya et al. (2002) was too poor (16 km s~^). Hence, we 
have re-analyzed the GS(3-2) data by Gesaroni et al. (2003), as the GS emission fine was 
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much stronger and observed with sufficient spectral resolution (0.5 km s~^). In this case the 
line profile deviates significantly from a Gaussian, presenting prominent emission in the red 
wing. Therefore we preferred to estimate the velocity from the first moment^ computed over 
a velocity interval including only the peak of the emission, from 108 to 116 km s~^, rather 
than from a Gaussian fit. The resulting Vlsr map is shown in Fig. 2a. 

The first conclusion that can be drawn from this figure is that all cores show clear 
velocity gradients, with Vlsr increasing steadily along well defined directions. We examine 
three possible explanations for such gradients: expansion, infall, or rotation. The first can be 
ruled out as the velocity gradient should be maximum in the same direction as the molecular 
outfiow, which is clearly not the case (see Figs. 2b and 2e). Spherical infall is also impossible, 
because self- absorption would shift the peak velocity towards lower values at the core center, 
whereas we observe a steady velocity increase along a well defined direction (see Figs. 2a, 
2c, 2d, and 2e). The fact that such a direction is perpendicular to the outflow axis strongly 
favours the rotation hypothesis: this is exactly what one expects if the core is rotating 
about the axis of the corresponding outfiow. This behaviour might be mimicked also by 
two distinct cores with different Vlsr and too close to be resolved by our observations: the 
velocity gradient would be a consequence of line emission from the two cores observed in 
the same instrumental beam. However, we believe this to be very unlikely. In fact, at least 
in the case of G31, the angular separation between the regions emitting at the maximum 
and minimum velocities is definitely greater than the beam size. As for G24, the two cores 
Al and A2 likely correspond to two distinct rotating disks. However, one cannot rule out 
another explanation, namely that they are part of the same geometrically thick rotating disk: 
in this case, the two emission peaks would be produced by the interaction of the outflow 
with the dense material of the disk. 

In conclusion, we believe that the most plausible explanation for the kinematics of G24 
and G31 is that the cores have toroidal structures undergoing rotation about the 
corresponding outflow cixis. Hereafter, we shall refer to these simply as disks, although 
one has to keep in mind that these are very different from the geometrically thin circumstellar 
disks seen in low-mass YSOs. 



°the first moment between vi and V2 is defined as }12t d ' where Tb is the brightness temperature. 
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3.3. Nature of the G24 and G31 disks 

The major question raised by our results is whether the disks are stable entities. In 
Table 1 we give a few disk parameters, among which the mass of the cores, Mgas, and the 
dynamical mass, Mjyn, needed for equilibrium. The former was estimated from the millimeter 
continuum emission assuming a mass opacity of ~ 0.02 cm~^ g~^ at 1.4 mm for a gas-to-dust 
ratio of 100 (see, e.g., Andre, Ward-Thompson, & Barsony 2000), and the temperatures 
listed in Table 1; the latter was computed assuming equilibrium between centrifugal and 
gravitational forces from the expression Afdyn = vf.^tRs\.v? i/G, where fj-ot is the rotation 
velocity, R is the radius of the disk, and i is the inclination angle of the disk assumed to be 
45°. Mgas is much larger than Mdyn, suggesting that the disks may be unstable. In principle, 
magnetic fields could stabilize the disks, but this would require a few 20-40 mG, values too 
large to be plausible even in regions as dense as 10^ cm~^ (see Fig. 1 of Crutcher 1999). 
Therefore, the disks must be transient structures with lifetimes of the order of the free-fall 
time, ts, also listed in Table 1. Another estimate of the disks hfetime, t^^, can be derived 
from the ratio between the disk mass, Mgas, ^-nd the accretion rate. The latter is computed 
from the expression Mace = 27rSit!i'in where E = Mgas/nR'^ is the surface density and t'in 
is the infall velocity, which has been assumed to be equal to the rotation velocity following 
Allen, Li, & Shu (2003). As one can see from Table 1, is very close to taco and both 
agree within a factor <4 with the outflow age, tout) derived from the data of Olmi et al. 
(1996b) and Furuya et al. (2002). Note that tout is to be multiplied by cot^^ to correct for 
the (unknown) inclination angle 6 of the flow with respect to the line of sight. We believe 
that 9 cannot differ signiflcantly from 45° otherwise blue- and red-shifted emission would 
mix up in the plane of the sky (for 9 ~ 90°) or along the fine of sight through the center (for 
9 ~ 0°): this implies a correction factor of order unity. The correction factors for 9 — 30° 
and 60° would be 1.7 and 0.6 respectively. 

In conclusion, the hfetime of the disks seems to be of order of 10^ yr. Such a short lifetime 
should imply 10 times less disks than UC Hll regions, which are supposed to live 10^ yr 
(Wood & Churchwell 1989). Although it is obviously impossible to confirm this estimate on 
a statistical ground, disks appear to be an ubiquitous phenomenon in massive star forming 
regions, as we have detected 4 of them in 2 regions only. Therefore, it seems unlikely that 
disks are 10 times less numerous than UC Hll regions. This implies a significantly longer 
lifetime than ~10^ yr, which in turn means that disks must be fed by a larger scale reservoir 
of material at a rate comparable to Mace — 10^^ Mq yr~^. Accretion rates that large have 
been estimated by Fontani et al. (2002) for the parsec-scale clumps where high-mass star 
formation is observed. 

If the disk lifetime is comparable to that of UC Hll regions, then the total accreted 
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mass should result in Mace x lOte ~ 10~^ x 10^ = 10^ M© of stars, too large a value for a 
single star, but acceptable if the infalling gas is accreting onto a cluster of stars. Indeed, this 
resembles the situation in the Orion cluster (Palla & Stabler 1999), where only <7% of the 
mass in stars (>600 Mq) is contained in the most massive star of the cluster (~40 Mq). 

The main finding of our study is that we have detected rotating disks associated with 
high-mass YSOs, and hence, this result strongly suggests that non-spherical accretion is a 
viable mechanism to form high-mass stars. Only a larger number of observations may confirm 
this conclusion on a statistical ground, thus assessing that disks are a natural product of the 
star formation process also for early-type stars. 

It is a pleasure to thank the staff of IRAM for his help during the observations. We also 
thank Daniele Galli for stimulating discussions about models of massive accretion disks. 
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Fig. 1. — Upper panel: Overlay of the 1.4 mm continuum contour map on the gray 
scale image obtained integrating the CH3CN(12-11) emission under the = 0, 1, and 2 
components for G24. Contour levels range from 0.02 to 0.2 in steps of 0.06 Jybeam^^. 
Grayscale levels range from 0.1 to 1.00 in steps of 0.18 Jy beam^^ km s^^. The straight 
line represents the axis of the bipolar outflow. The synthesized beam is shown in the lower 
right-hand corner. Lower panel: same as upper panel for G31. Gontour levels range from 
0.08 to 1.28 in steps of 0.3 Jybeam"^. Grayscale levels range from 0.1 to 0.82 in steps of 
0.12 Jybeam~^km s~^. 
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Fig. 2. — a. Map of the first moment of the CS(3-2) hne observed by Cesaroni et al. 
(2003) towards G24 C. The contour levels in km s~^ are indicated in the wedge to the right 
of the panel. The position offsets are relative to a(J2000) = 18^^36°^12!66, 5(J2000) = 
— 07°12'10'.'15. The straight line represents the outflow axis. The star indicates the peak 
of the mm continuum emission, b. Comparison between the bipolar outflows observed 
by Furuya et al. (2002) {contours) , the CH3CN(12-11) line emission mapped by us towards 
G24 Al and G24 A2 {gray scale), and the CS(3-2) line map from Furuya et al. (2002) towards 
G24 C {gray scale). Grayscale levels for CH3CN(12-11) are the same as in Fig. 1, while for 
CS(3-2) grayscale levels range from 0.4 to 0.8 in steps of 0.1 Jybeam~^km s~^. Blue and 
red contours represent respectively the blue- and red-shifted ^^CO(l-O) line emission (see 
Fig. 2 of Furuya et al. 2002). c. Map of the CH3CN(12-11) line peak velocity towards 
G24 A2 obtained with a Gaussian fit. The contour levels in km s~^ are indicated in the 
wedge to the right of the panel. The position offsets are relative to q;(J2000) = 18'^36™12!66, 
5(J2000) = — 07°12'10'.'15. The dashed line indicates the direction of the outflow, arbitrarily 
associated with G24 Al (see text). The star indicates the peak of the mm continuum 
emission, d. Same as c for G24 Al. e. Same as c for G31. The position offsets are 
relative to a(J2000) = 18^47°'34^33, (5(J2000) = -01°12'46'.'50. Note that for this case the 
direction of the outflow has been obtained connecting the peaks of the blue and red lobes, 
f. Comparison between the bipolar outflow observed by Olmi et al. (1996b) {contours), 
and the CH3CN(12-11) line emission map towards G31 {gray scale). Grayscale levels for 
CH3CN(12-11) are the same as in Fig. 1. Blue and red contours represent respectively the 
blue- and red-shifted CO (1-0) line emission (see Fig. 5 of Olmi et al. 1996b). 
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Table 1. Parameters of disks and outflows in G24 and G31. 
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G24 Al 


80*' 


0.02 


1.50 


23 
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4 X 10^ 


7 X 10^ 


2 X 10"' 


2 X 10" 


-2 


5 X 10" 


-4 c 


G24 A2 


so'' 


0.02 


0.75 


4 


80 


3 X 10^ 


1 X 10"' 




8 X 10" 


-3 






G24 C 


so'' 


0.04 


0.50 


5 


250 


9 X 10^ 


4 X lO"* 


2 X 10" 


6 X 10" 


-3 


5 X 10" 


-4 c 


G31 


230'* 


0.04 


2.10 


87 


490 


6 X 10^ 


1 X 10* 


2 X 10^ 


5 X 10" 


-2 


4 X 10" 


-4 d 



^Computed assuming Vin=VTot- 
^ Prom Codella et al. (1997). 
" Prom Puruya et al. (2002). 
d Prom Olmi et al. (1996b). 



